Direct observation of large temperature fluctuations during DNA thermal 
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In this paper we report direct measurement of large low frequency temperature fluctuations in 
double stranded (ds) DNA when it undergoes thermal denaturation transition. The fluctuation, 
which occurs only in the temperature range where the denaturation occurs, is several orders more 
than the expected equilibrium fluctuation. It is absent in single stranded (ss) DNA of the same 
sequence. The fluctuation at a given temperature also depends on the wait time and vanishes in a 
scale of few hours. It is suggested that the large fluctuation occurs due to coexisting denaturatcd 
and closed base pairs that are in dynamic equilibrium due to transition through a potential barrier 
in the scale of 25 - 30k B T (T = 300K). 
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INTRODUCTION 



The physics of thermal denaturation has been a topic 
of investigation for the past four decades 0, 0, S 0] • 
In thermal denaturation, which we call melting transition 
(also known as helix-coils transition) the hydrogen bonds 
between the base pairs of the double stranded DNA (ds- 
DNA) break and the dsDNA separates into two single 
stranded DNA (ssDNA). It occurs over a range of tem- 
perature and there are distinct steps which depend on 
the sequence as well as the length of the DNA (number 
of base pairs) [|| . The AT segments with lower binding 
energy melt at lower temperatures compared to the GC 
segments. The melting of DNA starts with creation of 
denaturation bubbles. As the denaturation process pro- 
ceeds the number of bound pairs decrease. 

In this paper we address the issue of temperature 
fluctuations during the process of thermal denaturation 
which we believe has not been addressed before. We 
present the first direct measurement of the temperature 
fluctuations in a number of double stranded DNA under- 
going denaturation. These include heterogeneous poly- 
mers like a 40 kbps T7 bacteriophage, 7250 bps plasmid 
(both ssDNA and dsDNA), lOObp heteropolymeric ds- 
DNA fragment and a 100bps homopolymer (...ATAT...). 
Our experimental results show that the process of ther- 
mal denaturation is associated with large low frequency 
fluctuations. These grow as the temperature enters the 
melting regime and die down once the melting process is 
complete. At a given temperature in the melting regime, 
the increased fluctuation shows a time dependence. It 
first grows and then decays. The basic process of ther- 
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FIG. 1: (a) Schematic of the calorimeter, (b) schematic of 
the noise measurement setup. 



mal detection is explained in the schematic of figure 1. 
The calorimeter consists of a substrate (heat capacity 
C su b ~ 30m J / K) and a thin film of Pt as a thermometer. 
The calorimeter is connected to a temperature controlled 
copper thermal bath (heat capacity Cbath ~ 60 J/ K) 
through a link of thermal resistance Rth ~ SOCff/V^.The 
total heat capacity of the system is the heat capacity 
of the sample C s plus that of the substrate C su b- But 
C s being much less (~ 10~ 6 J/K), the total heat capac- 
ity is ~ C' su b. The calorimeter has a relaxation time 
t = RthCsub- We use two calorimeters connected in 
a bridge arrangement which enhances the sensitivity of 
the experiment and also reduces the effect of small ther- 
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FIG. 2: (a) Cartoon of the three contributions to the 
observed Sv(f)- (b) The observed Sv(f) and S t h at 
T — 307K(« temperature where the denaturation sets in) 
for empty calorimeter (open symbols) and calorimeter with 
7250bp dsDNA (closed symbols) (c) the observed Sv{f) and 
Sth at T — 326K for empty calorimeter (open symbols) and 
calorimeter with 7250bp dsDNA (closed symbols). 



mal drifts. The sample (DNA) is mounted on the sub- 
strate. When the strands separate out the DNA absorbs 
energy from the calorimeter leading to lowering of the 
calorimeter temperature in a time scale that is faster 
than the thermal relaxation time (r) of the calorime- 
ter. At a longer timescale the calorimeter equilibrates 
with the heat bath. In the timescale < r, if the energy 
exchanged by the DNA with the calorimeter/substrate 
is AU, the temperature change AT of the calorimeter 
is AT w AU/C su h- As the denaturation transition pro- 
ceeds, the process of energy exchange thus leads to a 
fluctuation in the temperature of the calorimeter (about 
its average temperature). The temperature change AT 
of the calorimeter arising from the energy exchange leads 
to a resistance change AR of the thermometer. The ther- 
mometer being current biased, this is detected as a volt- 
age change (AV). The voltage fluctuation thus arising 
from the temperature fluctuation is given by: 



{(AV) 2 ) 
V 2 



(J3T) 



2 ((AT) 2 ) 

rp2 



(1) 



where V is the voltage across the Pt thermometer. (3(= 
J^Tj^) of the platinum thermometer was determined from 
the calibration of the thermometers and « 3.8 x 10~ 3 /K 
around 300K. The thermal link of the calorimeter deter- 
mines the bandwidth of the detection system which is 
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FIG. 3: An example of the PSD of the temperature fluctua- 
tion S T (Sam P ie) for 7250 bps dsDNA at T=326K (shown by 
cirles). The PSD of temperature fluctuation due to the tem- 
perature fluctuation of the bath is shown by a dotted line. 



like a low pass filter with a roll off at f c = jp— . 

The temperature fluctuations of the calorimeter can 
be characterized by the power spectral density (PSD) 
and the mean square temperature fluctuation ((AT) 2 ). 
The observed temperature fluctuation though finite is 
small and a digital signal processing technique Q was 
used to eliminate the different background contributions 
(described later on). Briefly, the voltage fluctuation 
(and from it the temperature fluctuation) was deter- 
mined by digitizing the AV of the thermometer as a 
function of time at a digitization rate of 1024 /sec. To 
reduce the contribution of the noise of the detector elec- 
tronics, the voltage fluctuation was measured using a 
frequency modulation technique where the fluctuation 
appears as a side band around a carrier frequency^ 
230Hz) 0- The data is digitally filtered and deci- 
mated and the PSD of voltage fluctuation Sv(f) was 
obtained by a fast fourier transform (FFT). As described 
later, from the observed Sy(f), the Nyquist noise and 
the 1// noise of the thermometer (predetermined exper- 
imentally using the empty calorimeter) were subtracted 
out to obtain the PSD of fluctuation from the sample 
Ssampie(f)- The mean square fluctuation ((AV~) 2 )/V 2 
= fn S S am P ie(f)df where f x and f 2 are the lower (w 
10~ 3 Hz) and upper (w 10Hz) frequency limits. The 
temperature fluctuation ((AT) 2 ) /T 2 was then obtained 
from ((AV) 2 )/V 2 by equation ^ The determination 
of ((AT) 2 ) /T 2 through PSD was necessary for proper 
subtraction of the contribution of Nyquist noise and the 
thermometer 1/f noise. 

A few /j,g to a few tens of /ig of a solution of the DNA in 
water are put on the substrate using a micropipette and 
the solution is allowed to evaporate off. DNA with sol- 
vent evaporated off harbours about 66% bound water Q 
and the conformation of the DNA does not undergo any 
change. We have observed similar fluctuations near the 
denaturation temperature when the experiment is carried 
with DNA in a buffer. We have used DNA with solvent 
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evaporated off because DNA in soiution reduces the sen- 
sitivity of the experiment due to the thermal mass of the 
solvent. 

Figure 2 shows the observed Sv(f) across the Pt film 
taken with an empty calorimeter as well as that taken 
with a sample of 7250bp dsDNA. It can be seen from 
figure 2(b) that at T = 307 K both the data are nearly the 
same. However at T — 326JT (figure 2(c)) the observed 
fluctuation with DNA is much larger. This is an example 
of a typical data. The observed Sv(f) can be written 
(shown in figure 2(a)) as: 



Sv(f) = S t h + Spt + S s 

ample 



(2) 



Here S t h is the white Nyquist noise given by Ak^TR 
(rs 10~ 1S V 2 / Hz) and S p t is predominantly the 1// noise 
of the Pt film resistor. The sum Sth + S p t (that con- 
stitutes the background) is measured by doing the runs 
on an empty calorimeter and is subtrated out from the 
measured Sy(f) to obtain S sam pie- (Note: Sth can be 
independently obtained by digitizing the out-of- phase 
signal from the lock- in amplifier |7J.) 

In figure 3 we show an example of the PSD of the tem- 
perature fluctuation St (Sample) (obtained from S samp i e 
using the relation S T ( sample ) (/) = ((iV)' 2 S samp i e (f)) for 
the 7250 bp dsDNA at T=326K. This figure shows that 
there is a large temperature fluctuation that has a pre- 
dominant low frequency component. In the same figure 
we also show the expected contribution from the bath 
fluctuation by a dotted line |!j . It can be seen that this 
is several orders smaller than the fluctuations due to the 
sample. 

The data were taken using two different procedures. 
In the first procedure, the temperature was stabilized at 
a particular value and two data sets were taken. Imme- 
diately the temperature was raised to the next desired 
value. At temperatures below that of the onset of melt- 
ing, the observed voltage fluctuation is mostly the 1/f 
noise of the thermometer. By slowly increasing the tem- 
perature we find the onset temperature where there are 
additional contributions to the fluctuations. Using this 
procedure we obtained data as a function of temperature 
for fixed wait time t w at every temperature. In the sec- 
ond method, the thermal bath was stabilized at a tem- 
perature where one sees significant fluctuations (which 
is known from data got using the first procedure). Data 
were taken at that temperature at different time intervals 
until the additional fluctuations fell below the detection 
limit. This way, using different wait time t w we can follow 
the approach to equilibrium. 

In figure 4 we show the ((AT) 2 ) /T 2 as a function of 
temperature for the different DNA samples. The figure 
shows that in a certain region in temperature the ds- 
DNA shows a large fluctuation which is absent in the ss- 
DNA with the same sequence. We find that ((AT) 2 )/T 2 
can be rather large, reaching close to lppm in the tem- 
perature range of interest . In the inset of figure 4 we 
show the UV absorption data for the T7 DNA. We can 
see that the melting temperature matches very well with 
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FIG. 4: The relative mean temperature fluctuation ^^ffl - 
as a function of temperature in the different DNA samples. 
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FIG. 5: (a) ^^ffl ^ as a function of the wait time t w at 329K 
for the T7 DNA. (b) the PDF of temperature fluctuation as 
a function of the wait time t w at 329K for the T7 DNA. The 
data are offset for clarity. The deviation from the Gaussian 
behavior is seen as deviation from the straight lines. 



that of increased fluctuation. Thus we can relate the 
observed temperature fluctuation in the dsDNA to the 
process of thermal denaturation. The fluctuation is ab- 
sent both in the fully zipped and fully unzipped state. 
The range in temperature where the excess noise persists 
depends on the length and heterogeneity, being maxi- 
mum for the T7 DNA (ST 15K ) and smallest for the 
100bps homopolymer (ST m 5K). It can also be seen that 
while in the lOObp heteropolymer the fluctuation jumps 
by more than an order of magnitude during melting, in 
the lOObp homopolymer the change is much lesser. In 
the case of the T7 DNA, the fluctuation during melting 
is about three orders higher than that of the background 
fluctuation. It is clear from these observations that both 
length and heterogeniety play important roles in deter- 
mining the magnitude of fluctuation. To summarize, the 
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fluctuation is larger in longer DNA samples and in het- 
erogenous samples. 

Figure 5(a) shows the '- as a function of the wait 

time t w at 329K for the T7 DNA. We define wait time 
as the time elapsed from when the temperature reached 
a stable value. This is a typical data set at the on- 
set of the fluctuation. Initially the fluctuation is very 
small. It increases with time, reaches a maximum and 
then decreases again. The t w dependence of the fluctu- 
ation is also clearly seen in the probability distribution 
function (PDF) of the temperature jumps (AT) shown 
in figure 5(b). The PDF is shown as relative probability 
vs (AT) 2 . The data are offset for clarity. For Gaussian 
distribution these plots should be straight lines. We see 
that at t w — the plot is a straight line. With larger t w 
there is deviation from straight line. This deviation in- 
creases with t Wl reaches a maximum and then decreases 
and finally the PDF again becomes Gaussian. For a given 
t w the data are shown both for the positive and negative 
AT. We note that for large AT a small yet definite asym- 
metry develops in the distribution function |T^ . 

During the denaturation process the observed tem- 
perature fluctuation ((AT) 2 ) /T 2 is much larger than 
what one would expect from the equilibrium fluctua- 
tion. For a microcanonical ensemble at temperature T 
and heat capacity C the temperature fluctuation is given 

by <( y! )2> = For the calorimeter at 300^ the esti- 
mated relative fluctuation « 4.5 x 1CP 22 . This is several 
orders less than the observed fluctuation. 

We would like to suggest that the observed large fluctu- 
ation is a result of coexisting phases that are in dynamical 
equilibrium. The transition between two phases has en- 
ergetics associated with it because it involves making and 
breaking of bonds leading to exchange of energy with the 
calorimeter. The observed slow dynamics suggests that 
the transition between the two phases involves thermal 
activation across barriers. The co-existing phases can 



arise from a number of sources like co-existence of loops 
(denaturated) and bound segments or due to the sys- 
tem being trapped in mismatched configuration. At this 
time we cannot definitely identify what leads to these co- 
existing phases observed in the experiment. We can get 
an estimate of the average activation energy AE from 
the wait time (t w ) dependence of the fluctuation using 
r = T exp(\AE\/k B T). Here r ~ 10 -7 s is the micro- 
scopic time for a basepair to move together or apart (in 
the absence of tension or base pairing interactions) |1C| . 
We find that AE ~ 25 - 30k B T where T = 300K. The 
value of AE is similar to the energy at the critical force 
of unzipping when the DNA is separated using force . 
We also note that our observation is distinct from what 
one expects from co-existing phases in a first order tran- 
sition. The observation of co-existing phases and the 
resulting fluctuation is metastable because it dies down 
after a long wait time. Here the system appears to be 
trapped for a long time with co-existing phases before it 
reaches thermal equilibrium. 

To summarize, we have made a direct measurement of 
temperature fluctuations during thermal denaturation of 
dsDNA. The fluctuation has a significant low frequency 
component and it peaks in the temperature region of 
denaturation. It is absent both in the fully closed and 
fully denaturated state. The fluctuation also has a wait 
time dependence and decays to zero after a sufficient wait 
time. We suggest that the fluctuation arises due to co- 
existing phases in dynamical equilibrium. 



II. ACKNOWLEDGMENT 

AKR and DC acknowledge financial support from DBT 
(Government of India) for a sponsored project and Na- 
gapriya acknowledges CSIR for support. 



[1] D. Poland, H. A. Scheraga, J. Chem. Phys. 45, 1456 

(1966) and Ibid 45, 1464 (1966). 
[2] R. M. Wartell and A. S. Benight, Phys. Rep. 126 67, 

(1985) 

[3] C. Richard and A. J. Guttmann, J. Stat. Phys. 115, 925 
(2004) 

[4] M. Barbi, S. Lepri, M. Peyrard, and N. Theodorakopou- 

los Phys. Rev. E 68, 061909 (2003) 
[5] T. Garel, C. Monthus and H. Orland, Euro. Phys. Letts. 

55 132, (2001) 

[6] Ralf Blossey and Enrico Carlon, Phys. Rev. E, 68, 061911 
(2003). 

[7] Arindam Ghosh, Swastik Kar, Aveek Bid and A. K. Ray- 
chaudhuri arXiv:cond-mat/0402130 vl 4 Feb 2004. J. 
H. Scoffield, Rev. Sci. Instr. 58, 985 (1987). 

[8] A.G.W. Leslie, S. Arnott, R. Chandrasekharan and RL. 
Ratliff, J. Mol. Biol. 143, 49 (1980) 

[9] The response of the calorimeter to the bath temper- 



ature fluctuations is that of a low pass thermal filter 
with roll off at /* = — ^ where r bs = R t hC ba th is 
the time constant of bath to DNA-substrate system. The 
contribution of the temperature fluctuation of the bath 
to the power spectral density (PSD) can be given by 



[10] S. Cocco, J. F. Marko and R. Monasson, Eur. Phys. J E 
10, 153 (2003) 

[11] Claudia Danilowicz, Vincent W. Coljee, Cedric 
Bouzigues, David K. Lubensky, David R. Nelson 
and Mara Prentiss, Proc. Natl. Sci. Acad., 100, 1694 
(2003) 

[12] The finite asymmetry in the PDF is a signature of non- 
equilibrium nature of the system. As the system moves 
towards equilibrium by absorbing heat from the sub- 
strate, there are more cooling jumps than heating jumps. 



